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FOREWORD 


Methodology  and  Characterization  Studies  during  fiscal  years 
1977  and  1978  included  gravity  surveys  in  ten  valleys  in  Arizona 
(five),  Nevada  (two).  New  Mexico  (two),  and  California  (one). 
The  gravity  data  were  obtained  for  the  purpose  of  estimating  al¬ 
luvial  gross  structure  and  shape  of  the  basins  and  the  thickness 
of  the  valley  fill.  There  was  also  the  possibility  of  detecting 
shallow  rock  in  areas  between  boring  locations.  Generalized 
interpretations  from  these  surveys  were  included  in  Fugro 
National's  Characterization  Reports  ( FN-TR- 25a  through  e)  . 

During  the  FY  77  surveys,  the  measurements  were  made  to  form  an 
approximate  one-mile  grid  over  the  study  areas,  and  contour  maps 
showing  interpreted  depth  to  bedrock  were  made.  In  FY  79,  the 
decision  was  made  to  concentrate  the  available  funds  on  the 
basic  Verification  Program  to  verify  and  refine  suitable  area 
boundaries.  This  decision  resulted  in  a  reduction  in  the 
gravity  program.  Instead  of  obtaining  gravity  data  on  a  grid, 
the  reduced  program  consisted  of  obtaining  gravity  measure¬ 
ments  along  profiles  across  the  valleys  where  Verification 
Studies  were  also  performed. 

The  Defense  Mapping  Agency  (DMA),  St.  Louis,  was  also  requested 
to  provide  gravity  data  from  their  library  to  supplement  the 
gravity  profiles.  For  Big  Smoky,  Reveille,  and  Railroad  val¬ 
leys,  a  sufficient  density  of  library  data  is  available  to 
permit  construction  of  interpreted  contour  maps  instead  of 
two-dimensional  cross  sections. 

In  late  summer  of  FY  79,  supplementary  funds  became  available  to 
begin  data  reduction.  At  this  time,  inner  zone  terrain  correc¬ 
tions  began  on  the  library  data  and  the  profiles  from  Big  Smoky 
Valley,  Nevada,  and  Butler  and  La  Posa  valleys,  Arizona.  The 
profile  data  from  Whirlwind,  Hamlin,  Snake  East,  White  River  and 
Garden  Coal  valleys,  Nevada  were  available  from  the  field  in 
early  October,  1979. 

A  continuation  of  gravity  interpretations  has  been  incorporated 
into  the  FY  80  contract  and  the  results  are  being  summarized  in 
a  series  of  valley  reports.  The  reports  covering  Nevada-Utah 
gravity  studies  will  be  numbered,  "FN-TR-33-",  followed  by  the 
abbreviation  for  the  subject  valley.  In  addition,  more  detailed 
reports  of  the  results  of  FY  77  surveys  in  Dry  Lake  and  Ralston 
valleys,  Nevada  are  being  prepared.  Verification  Studies  are 
continuing  in  FY  80  and  gravity  studies  are  included  in  the 
program.  DMA  will  continue  to  obtain  the  field  measurements  and 
it  is  planned  to  return  to  the  grid  pattern.  The  interpretation 
of  the  grid  data  will  allow  the  production  of  contour  maps  which 
will  be  valuable  in  the  deep  basin  structural  analysis  needed 
for  computer  modeling  in  the  Water  Resources  Program.  The 
gravity  interpretations  will  also  be  useful  in  the  Nuclear 
Hardness  and  Survivability  (NH&S)  evaluations. 
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The  basic  decisions  governing  the  gravity  program  are  made  by 
BMO  following  consultation  with  TRW  Inc.,  Fugro  National  and 
the  (DMA).  Conduct  of  the  gravity  studies  is  a  joint  effort 
between  DMA  and  Fugro  National.  The  field  work,  including 
planning,  logistics,  surveying,  and  meter  operation  is  done  by 
the  Defense  Mapping  Agency  Hydrographic/Topographic  Center 
(DMAHTC) ,  headquartered  in  Cheyenne,  Wyoming.  DMAHTC  reduces 
the  data  to  Simple  Bouguer  Anomaly  (see  Section  A1.4,  Appendix 
A1.0).  The  Defense  Mapping  Agency  Aerospace  Center  (DMAAC) ,  St. 
Louis,  calculates  outer  zone  terrain  corrections. 

Fugro  National  provides  DMA  with  schedules  showing  the  valleys 
with  the  highest  priorities.  Fugro  National  also  recommended 
locations  for  the  profiles  in  the  FY  79  studies  within  the 
constraints  that  they  should  follow  existing  roads  or  trails. 
Any  required  inner  zone  terrain  corrections  are  calculated  by 
Fugro  National  prior  to  making  geologic  interpretations. 
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1.0  INTRODUCTION 


1.1  OBJECTIVE 

Gravity  measurements  were  made  in  Hamlin  Valley  for  the  purpose 
of  estimating  the  overall  shape  of  the  structural  basin  and  the 
thickness  of  alluvial  fill.  The  estimates  will  be  useful  in 
modeling  the  dynamic  response  of  ground  motion  in  the  basin  and 
in  evaluating  ground-water  resources. 

1.2  LOCATION 

Hamlin  Valley  is  located  principally  in  Lincoln  and  White  Pine 
counties  in  east-central  Nevada  with  portions  in  Millard, 
Beaver,  and  Iron  counties,  Utah  (Figure  1).  Hamlin  Valley  is 
approximately  180  miles  (300  km)  NNE  of  Las  Vegas,  Nevada. 
U.S.  Highway  50  crosses  the  north  end  of  the  valley.  There  are 
no  paved  roads  within  the  valley,  but  access  is  generally  good 
along  well  maintained  ranch  and  mine  roads.  The  nearest  towns 
are  Baker,  Nevada  and  Garrison,  Utah,  along  Highway  73  near  the 
northwest  corner  of  the  valley. 

Hamlin  Valley  is  bounded  by  mountain  ranges  on  three  sides  and 
open  to  Snake  Valley  on  the  north  (Figure  2).  The  Burbank  Hills 


and 

the  Needle  Range 

fo  rm 

the  eastern  boundary  of 
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valley. 
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Rock 
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Most  of  the  valley  is  undeveloped  rangeland  with  several  ranches 
and  agricultural  fields  along  Hamlin  Wash. 
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1.3  SCOPE  OF  STUDY 

The  Defense  Mapping  Agency  Hyd rog raph i c/Topog raph i c  Center/- 
Geodetic  Survey  Squadron  (DMAHTO/GSS)  obtained  255  gravitational 
field  measurements  along  nine  cross  valley  profiles  in  Hamlin 
Valley.  Profile  positions  are  shown  in  Figure  2  and  the  loca¬ 
tions  of  the  individual  stations  are  shown  in  Figure  3.  The 
profiles  are  oriented  approximately  east-west.  They  are  between 
6  and  14  miles  (11  to  23  km)  long  and  are  spaced  approximately 
5  miles.  (8.5  km)  apart.  The  gravity  sampling  interval  was 
1  mile  (1.6  km)  over  the  central  valley  section  and  .25  mile 
(.4  km)  near  the  valley  boundaries.  The  more  dense  sampling  was 
used  on  the  valley  flanks  to  define  any  steep  gravity  gradients 
associated  with  boundary  faults  and  to  resolve  anomalies  with 
high  spatial  frequency  that  could  be  associated  with  shallow 
bed  rock . 


The  tolerance  for  establishing  the  station  elevations  was  5  feet 
(1.5  m)  ,  which  limits  the  gravity  precision  to  .3  milligals. 
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2.0  GRAVITY  DATA  REDUCTION 

DMAHTC/GSS  obtained  the  basic  observations  and  reduced  them  to 
Simple  Bouguer  Anomalies  (SBA)  for  each  station  as  described  in 
Appendix  A1.0.  Up  to  three  levels  of  terrain  corrections  were 
applied  to  convert  the  SBA  to  the  Complete  Bouguer  Anomaly 
(C BA).  First,  the  Defense  Mapping  Agency  Aerospace  Center 
(DMAAC) ,  St.  Louis,  used  its  library  of  digitized  terrain  data 
and  a  computer  program  to  calculate  terrain  corrections  out  to 
104  miles  (167  km)  from  each  station.  When  the  program  could 
not  calculate  the  terrain  effects  near  a  station,  a  ring  tem¬ 
plate  was  used  to  estimate  the  effect  of  terrain  within  approxi¬ 
mately  3000  feet  of  the  station.  The  third  level  of  terrain 
corrections  was  applied  to  those  stations  where  10  feet  or  more 
of  relief  was  observed  within  130  feet.  In  these  cases,  the 
elevation  differences  were  measured  in  the  field  at  a  distance 
of  130  feet  along  six  directions  from  the  stations.  These  data 
were  used  to  calculate  the  effect  of  the  very  near  relief.  The 
CBA  data  for  the  Hamlin  Valley  stations  are  listed  in  Appen¬ 
dix  A2.0. 
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3.0  GEOLOGIC  SUMMARY 

The  structural  geologic  setting,  major  rock  types,  and  deposi- 
tional  regime  of  the  valley  fill  material  are  important  consid¬ 
erations  in  the  interpretation  of  gravitational  field  data. 

Hamlin  Valley  is  an  elongate  north-south  trending  alluvial  basin 
exhibiting  typical  Basin  and  Range  structure.  Block  faulting 
has  formed  this  strong  north-south  physiographic  framework. 
According  to  Hintze  (1963),  a  fault  concealed  by  alluvium 
parallels  the  eastern  margin  of  the  valley.  Rocks  in  the  Needle 
Range  are  cut  by  predominantly  north-south  trending  faults. 
Most  faults  in  the  Snake  Range,  White  Rock  Mountains,  and 
Limestone  Hills  have  northwest  or  north-south  orientations. 

There  is  evidence  of  recent  upwarping  of  the  southern  end  of  the 
valley.  The  southwestern  lake  beds  are  at  a  relatively  high 
topographic  level  and  are  being  actively  eroded,  and  the  drain¬ 
age  trends  anomalously  to  the  north.  The  alluvium  is  offset  by 
numerous  small  faults  in  this  part  of  the  valley  (Fugro  Nation¬ 
al  ,  1979)  . 

The  mountains  in  the  northern  section  of  the  valley  are  composed 
of  an  assemblage  of  Paleozoic  limestones,  dolomites,  shales,  and 
sandstones  (Hintze,  1963;  Hose  and  Blake,  1976;  Tschanz  and 
Pampeyan,  1970).  The  White  Rock  Range,  on  the  south  and  the 
central  portion  of  the  Needle  Range,  on  the  east,  are  composed 
of  undifferentiated  Tertiary  volcanic  rocks. 
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Most  of  the  surficial  deposits  in  Hamlin  Valley  are  late  Terti¬ 
ary  and  Quaternary  alluvial  fan  and  lake  deposits.  Basin  fill 
deposits  are  described  in  the  Verification  Studies  ( FY  79, 
FN-TR-27-1A) .  Three  relative  ages  of  alluvial  fans  are  recog¬ 
nized.  The  older  fan  deposits,  consisting  of  well-indurated 
gravels  and  sandy  gravels  near  the  Needle  and  Snake  Ranges,  are 
the  least  extensive  unit.  Intermediate  age  alluvial  fan  de¬ 
posits  are  the  most  extensive  surficial  unit.  They  are  composed 
chiefly  of  silty  to  clayey  sands.  Younger  alluvial  fan,  ter¬ 
race,  and  modern  stream  channel  deposits  consist  dominantly  of 
silty  sand,  and  generally  occcur  in  the  central  part  of  the 
valley.  Lacustrine  deposits,  primarily  clayey  sands,  underlie 
about  one-fifth  of  the  valley  and  are  generally  restricted  to 
the  west  side  of  Hamlin  Wash  (Fugro  National,  1979). 
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4 . 0  INTERPRETATION 

A  negative  gravity  anomaly  over  a  valley  such  as  Hamlin  Valley 
is  created  when  low  density  alluvial  material  filling  the  valley 
is  surrounded  by  dense  rock  in  the  mountains.  Gravity  profiles 
across  such  a  valley  are  often  U-shaped,  low  in  the  middle  where 
the  fill  is  thickest,  and  high  on  the  ends  where  the  fill  thins 
and  disappears.  Interpretation  requires  removal  of  regional 
trends  leaving  the  gravity  reflection  of  the  valley  fill.  The 
gravity  data  and  interpreted  geologic  model  for  the  nine  pro¬ 
files  across  Hamlin  Valley  are  shown  in  Figures  4  through  12. 

4 . 1  REGIONAL  ~  RESIDUAL  SEPARATION 

A  fundamental  step  in  gravity  interpretation  is  isolation  of  the 
portion  of  the  C BA  which  represents  the  geologic  feature  of 
interest,  in  this  case  the  relatively  low  density  valley  fill. 
The  portion  of  the  CBA  which  corresponds  to  this  alluvial 
material  is  called  the  'residual  anomaly'. 

The  residual  anomaly  was  isolated  by  first  estimating  the  way 
the  CBA  field  would  have  appeared  if  there  had  been  no  valley 
fill  present.  This  estimated  field  is  called  the  'regional' 
gravity.  The  regional  gravity  is  subtracted  from  the  CBA  to 
produce  the  residual  anomaly.  For  this  study,  the  regional 
field  was  calculated  by  linear  interpolation  between  the  CB* 
values  at  bedrock  stations  on  opposite  ends  of  the  profiles. 
Where  only  one  end  of  a  profile  was  on  bedrock,  the  regional 
value  on  the  other  end  was  assigned  a  quantity  consistent  with 
the  regional  trend  in  the  area.  This  separation  technique  is 
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only  approximate.  Some  regional  effects  may  still  remain  after 
the  subtraction  but  the  error  is  probably  small  compared  to  the 
large  residual  anomaly  values  on  these  profiles. 

The  regional  field  used  for  each  profile  is  shown  together  with 
the  CBA  in  the  top  portion  of  Figures  4  through  12.  The  resid¬ 
ual  values  along  each  profile  are  shown  by  the  crosses  (x)  in 
the  center  portions  of  Figures  4  through  12. 

4.2  DENSITY  SELECTION 

The  construction  of  a  geologic  model  to  account  for  the  residual 
anomaly  requires  selection  of  density  values  representative  of 
the  basin  fill  and  of  the  underlying  rock.  Since  only  very 
generalized  density  information  is  available,  the  geologic 
interpretation  of  the  gravity  data  can  only  be  a  coarse  approxi¬ 
mation.  Average  in  situ  density  of  the  alluvial  fill  material 
was  measured  between  a  depth  of  100  and  160  feet  by  six  shallow 
borings  in  Hamlin  Valley.  The  observed  density  range  for  the 
soil  was  2.0  to  2.3  g/cm^.  These  borings  were  drilled  during 
the  Verification  Studies  of  Hamlin  Valley  (FY  79,  FN-TR-27-IV) . 
The  larger  density  value  was  used  in  the  modeling  process  to 
approximate  the  overall  density  increase  in  the  alluvium 
due  to  compaction  with  depth  (see  Grant  and  West,  1965). 

Published  values  for  carbonate  rocks  typically  range  between 
2.6  and  2.8  g/cm^.  The  Paleozoic  carbonate  rocks  in  Nevada 
are  generally  reported  to  be  relatively  high  in  density,  on  the 
order  of  2.8  g/cm^.  The  Nevada  volcanic  rocks  are  highly 
variable  in  density,  ranging  between  2.2  and  2.5  g/cm^. 
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The  calculated  basin  depths  are  very  dependent  on  the  density 
values  assigned  to  the  various  valley  materials.  A  one  percent 
change  in  the  average  alluvial  fill  density  will  result  in  a 
five  percent  change  in  the  calculated  fill  thickness. 

4 . 3  MODELING 

An  iterative  computer  program  that  calculates  the  gravitational 
field  for  two-dimensional  models  was  used  to  establish  a  thick¬ 
ness  of  alluvium  under  each  profile.  The  cross-sectional 
models  appear  as  a  set  of  either  1  or  0.5-km-wide  blocks 
whose  tops  are  at  surface  elevation  and  whose  bottoms  represent 
the  alluvium-bedrock  boundary.  The  elevations  at  the  bottoms  of 
the  blocks  were  adjusted  by  iterative  computation  until  the 
computed  gravity  anomaly  for  the  valley  fill  differed  by  less 
than  1  milligal  from  the  observed  residual  gravity  anomaly. 
The  computed  gravity  anomaly  from  the  final  model  is  shown  as  a 
continuous  line  in  the  center  portion  of  Figures  4  through  12. 
The  resulting  basin  models  are  shown  in  the  lowest  section  of 
Figures  4  through  12.  The  cross  sections  have  a  five  times 
vertical  exaggeration  so  that  gentle  slopes  appear  steep. 

The  gravity  survey  of  Hamlin  Valley  indicates  a  structural 
basin  which  was  formed  as  a  deep  graben  bounded  by  steep  faults. 
The  basin  appears  to  be  several  thousand  feet  deeper  in  the 
center,  near  profiles  HV-5  and  HV-6,  than  it  is  at  either 
end.  Actually,  it  is  effectively  terminated  near  profile 
HV-8.  Another  basin  appears  to  lie  southeast  of  profile  HV-8. 
Even  though  the  basin  is  elongate  N-S,  there  appears  to  be  a 
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component  of  E-W  deformation.  The  rock  outcrop  lines  strike 
approximately  E-W  for  about  6  miles  on  both  sides  of  the  valley. 
On  ’"he  east  side,  this  trend  occurs  near  profile  HV-2  and  on  the 
west  side  it  occurs  near  HV-3.  South  of  profile  HV-6,  the 
strike  of  the  rock  outcrops  is  NW-SE.  The  gravity  profiles 
indicate  that  the  axis  of  the  basin  also  shifts  in  direction,  or 
is  offset  laterally,  at  several  places  in  the  valley.  The  axis 
at  profile  HV-1  (which  may  actually  be  the  axis  of  Snake  Valley) 
appears  to  be  approx imately  5  miles  east  of  the  axis  at  profile 
HV-2,  and  it  apparently  trends  N-S  between  profiles  HV-2  and 
HV-3.  The  basin  appears  to  be  wider  and  more  complicated 
beneath  profiles  HV-4  and  HV-5  and  the  location  of  the  axis  is 
subject  to  question.  However,  it  seems  to  be  a  mile  or  so  west 
of  the  axis  at  HV-3.  The  orientation  of  profiles  HV-7  and  HV-8 
is  such  that  they  probably  do  not  show  the  axis  of  the  struc¬ 
tural  basin. 

Steep  gravity  gradients  in  this  basin  and  range  valley  are 
interpreted  as  being  caused  by  bedrock  faults.  The  faults  are 
interpreted  on  the  gravity  profiles  where  the  gradients  are 
maximum.  See  Figure  13  for  an  interpretation  of  possible  fault 
relationships  between  the  profiles.  Major-range  bounding  faults 
can  be  interpreted  with  confidence  on  profiles  HV-2  through 
HV-6 .  The  displacement  associated  with  the  boundary  faults 
increases  progressively  from  profiles  HV-2  through  HV-5,  and 
decreases  from  HV-5  to  HV-6.  If  these  faults  intersect  profiles 
HV-7  and  8,  they  cross  at  oblique  angles  and  show  only  a  compo¬ 
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Neither  of  the  boundary  faults  show  clearly  on  profile  HV-1. 
According  to  the  trend  established  from  profile  HV-5  to  HV-2, 
the  offset  at  HV-1  would  be  relatively  small  and  the  gravity 
expression  would  be  correspondingly  subtle.  If  the  suggestion 
that  there  has  been  E-W  deformation  between  profiles  HV-1  and 
HV-2  is  correct,  HV-1  may  have  crossed  only  the  west-boundary 
fault. 

Hintze  (1963)  shows  a  fault  covered  by  alluvium  along  the 
eastern  boundary  of  the  valley.  However,  the  gravity  inter¬ 
pretation  indicates  that  there  is  a  pediment-like  feature  on  the 
east  side  of  the  valley  and  that  rock  extends  westward  from  the 
outcrop  at  shallow  depths  for  2  to  4  miles  (3  to  7  km)  before 
being  faulted  downward.  No  corresponding  feature  is  shown  on 
the  western  valley  flanks,  where  the  boundary  fault  appears  to 
be  much  nearer  the  rock  outcrop  line. 

Profiles  HV-4  and  HV-5  indicate  that  there  may  be  one  or  more 
smaller  faults  within  the  basin,  forming  a  local  grAben  along 
the  western  boundary.  The  interpretation.  Figure  13,  shows  a 
fault  roughly  parallel  to  and  approximately  1  mile  (1.6  km) 
east  of  the  western  boundary  fault. 

The  gravity  survey  shows  the  southeastern  end  of  the  valley  to 
suddenly  become  substantially  shallower  at  Profile  lines  7  and 
8.  The  maximum  depth  of  the  basin  shown  on  these  profiles  varies 
from  about  1300  to  500  feet  (397  to  153  m)  .  The  alluvial  fill 
is  very  thin  at  the  center  of  each  of  these  profiles  suggesting 
an  E-W  trending  bedrock  ridge  extending  toward  the  extrusive 
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volcanic  rock  outcrops  mapped  in  this  area.  This  bedrock  ridge 
could  restrict  north-south  movement  of  ground  water.  An  alter¬ 
native  interpretation  of  profiles  HV-7  and  8  is  that  they  cross 
a  shallow  flow  of  igneous  material  with  a  density  greater  than 
the  alluvium.  If  this  is  the  case  there  may  be  substantial 
thickness  of  alluvium  beneath  the  flow. 

Profile  HV-9  crosses  the  southern  end  of  the  valley  where  it 
trends  NW-SE.  The  profile  extends  across  the  Needle  Range  into 
Pine  Valley.  The  pediment  feature  is  absent  here,  and  the 
gravity  interpretation  is  consistent  with  the  fault  suggested  by 
Hintze  (1963)  in  this  part  of  the  valley.  This  portion  of  the 
Needle  Range  is  interpreted  to  be  a  horst  since  a  fault  is 
interpreted  also  on  the  western  flank  of  Pine  valley.  A  bound¬ 
ary  fault  is  also  suggested  near  the  southwestern  end  of  profile 
HV-9.  It  is  assumed  that  this  fault  runs  parallel  to  the 
northern  boundary  of  the  White  Rock  mountains.  It  may  also  be 
seen  on  the  southern  end  of  profile  HV-8 .  The  depth  of  the 
basin  beneath  profile  HV-9  is  about  2800  feet  (853  m) . 
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5.0  CONCLUSIONS 

The  interpretation  of  the  gravity  survey  of  Hamlin  Valley 
indicates  that  there  are  major  range  bounding  normal  faults  on 
both  sides  of  the  valley.  The  major  graben  block  between  these 
boundary  faults  is  oriented  NNE-SSW.  It  is  calculated  to  be 
between  2450  feet  (747  m)  deep  in  the  north  end  of  the  valley 
and  6700  feet  (2042  m)  deep  near  the  valley  center. 

There  may  be  a  smaller  graben  trending  NNW  along  the  western 
side  of  the  major  block,  near  profiles  HV-4  and  HV-5.  This 
smaller  block  appears  to  be  from  1000  to  2000  feet  (305  to 
610  m)  deeper  than  the  main  block. 

A  bedrock  ridge,  which  could  restrict  all  but  the  shallowest 
ground  water  movement  is  interpreted  near  profiles  HV-7  and 
HV-8. 

There  is  a  large,  well  defined  negative  gravity  anomaly  as¬ 
sociated  with  Hamlin  Valley.  An  average  density  contrast  of 
0.50  g/cm^  between  the  alluvium  and  bedrock  was  used  to  cal¬ 
culate  the  thickness  of  the  valley  fill  material.  The  calcu¬ 
lated  bedrock  depth  can  only  be  an  approximation  since  little  is 
known  about  the  actual  density  distribution  in  and  around  the 
valley.  Future  studies  that  acquire  better  density  data  or 
depth  to  bedrock  in  deep  parts  of  the  valley  can  be  used  to 
refine  the  gravity  interpretation. 
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A1.0  GENERAL  PRINCIPLES  OF  THE  GRAVITY 
EXPLORATION  METHOD 


Al.l  GENERAL 

A  gravity  survey  involves  measurement  of  differences  in  the 
gravitational  field  between  various  points  on  the  earth's 
surface.  The  gravi tational  field  values  being  measured  are  the 
same  as  those  influencing  all  objects  on  the  surface  of  the 
earth.  They  are  generally  associated  with  the  force  which 
causes  a  1  gm  mass  to  be  accelerated  at  980  cm/sec^.  This 
force  is  normally  referred  to  as  a  1  g  force. 

Even  though  in  many  applications  the  g ravi tat ional  field  at  the 
earth's  surface  is  assumed  to  be  constant,  small  but  distin¬ 
guishable  differences  in  gravity  occur  from  point  to  point.  In 
a  gravity  survey,  the  variations  are  measured  in  terms  of  milli- 
gals.  A  milligal  is  equal  to  0.001  cm/second2  or  0.00000102  g. 
The  differences  in  gravity  are  caused  by  geometrical  effects, 
such  as  differences  in  elevation  and  latitude,  and  by  lateral 
variations  in  density  within  the  earth.  The  lateral  density 
variations  are  a  result  of  changes  in  geologic  conditions.  For 
measurements  at  the  surface  of  the  earth,  the  largest  factor 
influencing  the  pull  of  gravity  is  the  density  of  all  materials 
between  the  center  of  the  earth  and  the  point  of  measurement. 

To  detect  changes  produced  by  differing  geological  conditions, 
it  is  necessary  to  detect  differences  in  the  gravi tational  field 
as  small  as  a  few  milligals.  To  recognize  changes  due  to 
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geological  conditions,  the  measurements  are  "corrected"  to  ac¬ 
count  for  changes  due  to  differences  in  elevation  and  latitude. 

Given  this  background,  the  basic  concept  of  the  gravitational 
exploration  method,  the  anomaly,  can  be  introduced.  If,  instead 
of  being  an  oblate  spheroid  characterized  by  complex  density 
variations,  the  earth  were  made  up  of  concentric,  homogeneous 
shells,  the  gravitational  field  would  be  the  same  at  all  points 
on  the  surface  of  the  earth.  The  complexities  in  the  earth's 
shape  and  material  distribution  are  the  reason  that  the  pull  of 
gravity  is  not  the  same  from  place  to  place.  A  difference  in 
gravity  between  two  points  which  is  not  caused  by  the  effects  of 
known  geometrical  differences,  such  as  in  elevation,  latitude, 
and  surrounding  terrain,  is  referred  to  as  an  "anomaly." 

An  anomaly  reflects  lateral  differences  in  material  densities. 
The  gravi tational  attraction  is  smaller  at  a  place  underlain  by 
relatively  low  density  material  than  it  is  at  a  place  underlain 
by  a  relatively  high  density  material.  The  term  "negative 
gravity  anomaly"  describes  a  situation  in  which  the  pull  of 
gravity  within  a  prescribed  area  is  small  compared  to  the  area 
surrounding  it.  Low-density  alluvial  deposits  in  basins  such  as 
those  in  the  Nevada-Utah  region  produce  negative  gravity  anoma¬ 
lies  in  relation  to  the  gravity  values  in  the  surrounding 
mountains  which  are  formed  by  more  dense  rocks. 

The  objective  of  gravity  exploration  is  to  deduce  the  variations 
in  geologic  conditions  that  produce  the  gravity  anomalies 
identified  during  a  gravity  survey. 
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Al.  2  INSTRUMENTS 

The  sensing  element  of  a  LaCoste  and  Romberg  gravimeter  is  a 
mass  suspended  by  a  zero-length  spring.  Deflections  of  the 
mass  from  a  null  position  are  proportional  to  changes  in  gravi¬ 
tational  attraction.  These  instruments  are  sealed  and  compen¬ 
sated  for  atmospheric  pressure  changes.  They  are  maintained  at 
a  constant  temperature  by  an  internal  heater  element  and  thermo¬ 
stat.  The  absolute  value  of  gravity  is  not  measured  directly  by 
a  gravimeter.  It  measures  relative  values  of  gravity  between 
one  point  and  the  next.  Gravitational  differences  as  small  as 
0.01  milligal  can  be  measured. 

Al. 3  FIELD  PROCEDURES 

The  gravimeter  readings  were  calibrated  in  terms  of  absolute 
gravity  by  taking  readings  twice  daily  at  nearby  U3GS  gravity 
base  stations.  Gravimeter  readings  fluctuate  because  of  small 
time-related  deviations  due  to  the  effect  of  earth  tides  and 
instrument  drift.  Field  readings  were  corrected  to  account  for 
these  deviations.  The  magnitude  of  the  tidal  correction  was 
calculated  using  an  equation  suggested  by  Goguel  (1954): 

C  =  P  +  Ncos  0  (cos  0  +  sin  0  )  +  Scos  0  (cos  0  -  sin  0) 
where  C  is  the  tidal  correction  factor,  P,  N,  and  S  are  tine- 
related  variables,  and  0  is  the  latitude  of  the  observation 
point.  Tables  giving  the  values  of  P,  N,  and  S  are  published 
annually  by  the  European  Association  of  Exploration  Geophysi¬ 


cists 
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The  meter  drift  correction  was  based  on  readings  taken  at  a 
designated  base  station  at  the  start  and  end  of  each  day.  Any 
difference  between  these  two  readings  after  they  were  corrected 
for  tidal  effects  was  considered  to  have  been  the  result  of 
instrumental  drift.  It  was  assumed  that  this  drift  occurred  at 
a  uniform  rate  between  the  two  readings.  Corrections  for  drift 
were  typically  only  a  few  hundredths  of  a  milligal.  Readings 
corrected  for  tidal  effects  and  instrumental  drift  represented 
the  observed  gravity  at  each  station.  The  observed  gravity 
values  represent  the  total  gravitational  pull  of  the  entire 
earth  at  the  measurement  stations. 

Al. 4  DATA  REDUCTION 

Several  corrections  or  reductions  are  made  to  the  observed 
gravity  to  isolate  the  portion  of  the  gravitational  pull  which 
is  due  to  the  crustal  and  near-surface  materials.  The  gravity 
remaining  after  these  reductions  is  called  the  "Bouguer 
Anomaly."  Bouguer  Anomaly  values  are  the  basis  for  geologic 
interpretation.  To  obtain  the  Bouguer  Anomaly,  the  observed 
gravity  is  adjusted  to  the  value  it  would  have  had  if  it  had 
been  measured  at  the  geoid,  a  theoretically  defined  surface 
which  approximates  the  surface  of  mean  sea  level.  The  dif¬ 
ference  between  the  "adjusted"  observed  gravity  and  the  gravity 
at  the  geoid  calculated  for  a  theoretically  homogeneous  earth  is 
the  Bouguer  Anomaly. 
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Four  separate  reductions,  to  account  for  four  geometrical 
effects,  are  made  to  the  observed  gravity  at  each  station  to 
arrive  at  its  Bouguer  Anomaly  value. 

a.  Free-Air  Effect;  Gravitational  attraction  varies  inversely 
as  the  square  of  the  distance  from  the  center  of  the  earth. 
Thus  corrections  must  be  applied  for  elevation.  Observed 
gravity  levels  are  corrected  for  elevation  using  the  normal 
vertical  gradient  of: 

FA  =  -0.09406  mg/ft  (-0.3086  milligals/meter) 
where  FA  is  the  free-air  effect  (the  rate  of  change  of  gravity 
with  distance  from  the  center  of  the  earth).  The  free-air 
correction  is  positive  in  sign  since  the  correction  is  opposite 
the  effect. 

b.  Bouguer  Effect:  Like  the  free-air  effect,  the  Bouguer 

effect  is  a  function  of  the  elevation  of  the  station,  but  it 
considers  the  influence  of  a  slab  of  earth  materials  between 
the  observation  point  on  the  surface  of  the  earth  and  the 
corresponding  point  on  the  geoid  (sea  level).  Normal  practice, 
which  is  to  assume  that  the  density  of  the  slab  is  2.67  grams 
per  cubic  centimeter  was  followed  in  these  studies.  The  Bouguer 
correction  (Bc),  which  is  opposite  in  sign  to  the  free-air 
correction,  was  defined  according  to  the  following  formula. 

Bc  =  0.01276  (2.67)  hf  (milligals  per  foot) 

Bc  =  0.04185  (2.67)  hm  (milligals  per  meter) 
where  hf  is  the  height  above  sea  level  in  feet  and  hm  is  the 
height  in  meters. 
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c.  Latitude  Effect:  Points  at  different  latitudes  will  have 
different  "gravities"  for  two  reasons.  The  earth  (and  the 
geoid)  is  spheroidal,  or  flattened  at  the  poles.  Since  points 
at  higher  latitudes  are  closer  to  the  center  of  the  earth  than 
points  near  the  equator,  the  gravity  at  the  higher  latitudes  is 
larger.  As  the  earth  spins,  the  centrifugal  acceleration 
causes  a  slight  decrease  in  gravity.  At  the  higher  latitudes 
where  the  radii  of  the  circles  of  latitude  are  smaller,  the 
centrifugal  acceleration  diminishes.  The  gravity  formula  for 
the  Geodetic  Reference  System,  1967,  gives  the  theoretical  value 
of  gravity  at  the  geoid  as  a  function  of  latitude.  It  is: 

g  =  978.0381  (1  +  0.0053204  sin2  jz5  -  0.0000058  sin22tf)  gals 
where  g  is  the  theoretical  acceleration  of  gravity  and  0  is 
the  latitude  in  degrees.  The  positive  term  accounts  for  the 
spheroidal  shape  of  the  earth.  The  negative  term  adjusts  for 
the  centrifugal  accelerat ion . 

The  previous  two  corrections  (free  air  and  Bouguer)  have  ad¬ 
justed  the  observed  gravity  to  the  value  it  would  have  had  at 
the  geoid  (sea  level).  The  theoretical  value  at  the  geoid  for 
the  latitude  of  the  station  is  then  subtracted  from  the  adjusted 
observed  gravity.  The  remainder  is  called  the  Simple  Bouguer 
Anomaly  (SBA) .  Most  of  this  gravity  represents  the  effect  of 
material  beneath  the  station,  but  part  of  it  may  be  due  to 
irregularities  in  terrain  (upper  part  of  the  Bouguer  slab)  away 


from  the  station. 
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d.  Terrain  Effect;  Topographic  relief  around  the  station  has 
a  negative  effect  on  the  gravitational  force  at  the  station.  A 
nearby  hill  has  upward  gravitational  pull  and  a  nearby  valley 
contributes  less  downward  attraction  than  a  nearby  material 
would  have.  Therefore,  the  corrections  are  always  positive. 
Corrections  are  made  to  the  SBA  when  the  terrain  effects  were 
0.1  milligal  or  larger.  Terrain  corrected  Bouguer  values  are 
called  the  Complete  Bouguer  Anomaly  (CBA) .  When  the  CBA  is 
obtained,  the  reduction  of  gravity  at  individual  measurement 
points  (stations)  is  complete. 

Ai.  5  INTERPRETATION 

The  first  step  in  interpretation  is  to  separate  the  portion  of 
the  CBA  that  might  be  caused  by  the  lightweight,  basin-fill 
material  overlying  the  heavier  bedrock  material  which  forms  the 
surrounding  mountains  and  presumably  the  basin  floor.  Since  the 
valley-fill  sediments  are  absent  at  the  stations  read  in  the 
mountains,  the  CBA  values  at  these  bedrock  stations  are  used  as 
the  basis  for  constructing  a  regional  field  over  the  valley.  A 
regional  field  is  an  estimation  of  the  values  the  CBA  would  have 
had  if  the  light  weight  sediments  (the  anomaly)  had  not  been 
there . 

The  difference  between  the  CBA  and  the  regional  field  is 
called  the  "residual"  field  or  residual  anomaly.  The  residual 
field  is  the  interpreter's  estimation  of  the  gravitational 
effect  of  the  geologic  anomaly.  The  zero  value  of  the  residual 
anomaly  is  not  exactly  at  the  rock  outcrop  line  but  at  some 
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distance  on  the  "rock"  side  of  the  contact.  The  reason  for  this 
is  found  in  the  explanation  of  the  terrain  effect.  There  is  a 
component  of  gravitational  attraction  from  material  which  is 
not  directly  beneath  a  point. 

If  the  "regional"  is  well  chosen,  the  magnitude  of  the  residual 
anomaly  is  a  function  of  the  thickness  of  the  anomalous  (fill) 
material  and  the  density  contrast.  The  density  contrast  is  the 
difference  in  density  between  the  alluvial  and  bedrock  material. 
If  this  contrast  were  known,  an  accurate  calculation  of  the 
thickness  could  be  made.  In  most  cases,  the  densities  are  not 
well  known  and  they  also  vary  within  the  study  area.  In  these 
cases,  it  is  necessary  to  use  typical  densities  for  materials 
similar  to  those  in  the  study  area. 

If  the  selected  average  density  contrast  is  smaller  than  the 
actual  density  contrast,  the  computed  depth  to  bedrock  will  be 
greater  than  the  actual  depth  and  vice-versa.  The  computed 
depth  is  inversely  proportional  to  the  density  contrast.  A  ten 
percent  error  in  density  contrast  produces  a  ten  percent  error 
in  computed  depth.  An  iterative  computer  program  is  used  to 
calculate  a  subsurface  model  which  will  yield  a  gravitational 
field  to  match  (approximately)  the  residual  gravity  anomaly. 
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